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Abstract: We consider Lovelock theories of gravity in the context of AdS/CFT. We show 
that, for these theories, causality violation on a black hole background can occur well in the 
interior of the geometry, thus posing more stringent constraints than were previously found 
in the literature. Also, we find that instabilities of the geometry can appear for certain 
parameter values at any point in the geometry, as well in the bulk as close to the horizon. 
These new sources of causality violation and instability should be related to CFT features 
that do not depend on the UV behavior. They solve a puzzle found previously concerning 
unphysical negative values for the shear viscosity that are not ruled out solely by causality 
restrictions. We find that, contrary to previous expectations, causality violation is not 
always related to positivity of energy. Furthermore, we compute the bound for the shear 
viscosity to entropy density ratio of supersymmetric conformal field theories from d = 4 
till d = 10 - i.e., up to quartic Lovelock theory -, and find that it behaves smoothly as a 
function of d. We propose an approximate formula that nicely fits these values and has a 
nice asymptotic behavior when d — > oo for any Lovelock gravity. We discuss in some detail 
the latter limit. We finally argue that it is possible to obtain increasingly lower values for 
rj/s by the inclusion of more Lovelock terms. 

Keywords: AdS/CFT. Causality. Positive energy. Lovelock gravity. Conformal collider 
physics. Higher curvature corrections. Black holes. Shear viscosity. Instabilities. 
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1. Introduction 

The AdS/CFT correspondence has provided a tool to study hydrodynamical aspects of 
quantum field theories at strong coupling. This was particularly timely due to the advent 
of experiments that prompted the exploration of QCD in a region of phase space where it 
displays such behavior. One of the most striking predictions of AdS/CFT had to do with 
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the shear viscosity to entropy density ratio, 77/s, of strongly coupled plasmas. Interestingly 
enough, it was found that there is a universal value for this ratio, 77/s = 1/47T, for theories 
whose gravity dual is governed by the Einstein-Hilbert action [1] , regardless of the matter 
content, the number of supersymmetries, the existence or not of a conformal symmetry, 
and even the spacetime dimensionality. On the other hand, all measured values for this 
ratio in any quantum relativistic system are above this value. This led to speculations that 
v/ s — 1/47T might be an exact statement in quantum relativistic systems, the so-called KSS 
bound conjecture [2]. Indeed, it is possible to provide a hand waving argument, relying 
on a quasiparticle description of the plasma, that links the KSS bound to the Heisenberg 
uncertainty principle. This argument is questionable, though, since there are convincing 
hints supporting a non-quasiparticle description of strongly interacting plasmas. 

The KSS bound conjecture has been thoroughly scrutinized for many years (see, for 
example, [3] for a recent review). It turned out to be the case, however, that when quantum 
corrections are included in the gravitational action, under the form of curvature squared 
terms, the 77/s ratio can be smaller than the KSS bound [4, 5]. In particular, there are 
string theory constructions where this is the case [5, 6]. On general grounds, the coefficient 
of the curvature squared terms must be small. However, for the particular combination 
given by the Gauss-Bonnet (GB) term, one may consider finite values of the coefficient, 
as the gravitational theory is then two-derivative. In this case the ratio is modified to 
77/s = l/47r(l — 4A), A being the appropriately normalized GB coupling, and the KSS 
bound is violated whenever A is positive [4, 5]. Even if a priori the addition of a GB 
term would lead to an arbitrary violation of the KSS bound, it turns out that causality 
constraints arise preventing the possibility of going to arbitrarily low values of 77/s [7]. In 
the case of 4d CFTs, for instance, this imposes the constraint A < 9/100, which reduces 
the minimum value of 77/s by a factor of 16/25. 1 

The study of a possible bound is interesting both from a theoretical standpoint as 
well as from a phenomenological one - it has been found experimentally that the quark- 
gluon plasma created in relativistic heavy ion collisions appears to have a very low shear 
viscosity to entropy density ratio [9] . A similar result was also found recently in a radically 
different context: that of strongly correlated ultracold atomic Fermi gasses in the so-called 
unitarity limit [10]. Both systems have measured values of 77/s compatible with 1/4-7T. On 
the other hand, from a theoretical standpoint, the causality constraints coming from the 
behavior of the geometry near the boundary, that rule the attainability of lower values of 
77/s, have a beautiful holographic dual in the CFT side: they arise from positivity of energy 
conditions [11, 12, 13]. The perfect matching of these two quite fundamental restrictions on 
a physically sensible theory constitutes a striking check of the AdS/CFT correspondence. 

Whether this is still valid for higher-dimensional CFTs became a natural question that 
was subsequently answered in a series of papers [14, 15, 16]. Summarizing, it turned out 
that GB theories lead to a violation of the KSS bound in any spacetime dimensionality, and 
causality constraints exactly match positivity of energy bounds on the CFT side. The fact 
that this matching is valid regardless of the dimensionality is puzzling and seems to provide 

1 It is interesting to mention at this point that the link between causality violation and the viscosity 
bound does not apply for thermal theories undergoing a low temperature phase transition [8]. 
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clues on possible non-stringy versions of AdS/CFT. Our knowledge of higher-dimensional 
CFTs is, however, too poor yet to push these arguments forward. In higher dimensions 
one has in general the choice of including other Lovelock terms in the action. These are 
higher order in curvature that still lead to second order equations of motion [17]. A first 
step in the analysis of these theories was done in [18, 19]. There it was found that causality 
and positivity of energy constraints still match perfectly, provided a conjectural relation 
between field theory parameters and Lovelock couplings was true (something that is proved 
in the present article). Still, something new happens here. Since there are now more free 
parameters available, it turns out that causality alone cannot prevent the r//s ratio from 
becoming arbitrarily small (or even negative). 

In this paper we analyze in detail the fate of the rj/s ratio in AdS/CFT restricted to 
the case of Lovelock theories. We start by describing, in section 2, the lagrangians, AdS 
vacua and black brane solutions of such theories. The latter are determined by solving 
a simple polynomial of the same order as the highest Lovelock term in the action. We 
then establish the holographic dictionary for such theories in section 3, by computing the 
central charge Ct that controls the stress tensor 2-point function, and the two independent 
quantities, t<i and £4, appearing in the 3-point function [11]. With these parameters we can 
jump into the CFT side and establish positivity of energy bounds tantamount to a series 
of constraints on the couplings of the gravitational lagrangian. 

We then consider graviton fluctuations about the black hole background for a general 
Lovelock theory. In the high-momentum limit, one can recast their equations of motion as 
Schrodinger-type equations with surprisingly simple effective potentials, which also have 
the interpretation of the effective speed of light of the graviton at various radii. Near 
the boundary, demanding causality of graviton propagation leads to bounds which exactly 
match those of positivity of energy. Close to the horizon, in turn, we show that unstable 
modes can appear if the parameters of the theory are not tuned properly. This leads 
to another set of bounds discussed in section 4, which for the cubic Lovelock theory in 
seven dimensions prevents the rj/s ratio from becoming arbitrarily small. In the case of 
higher order Lovelock theories, a new phenomenon occurs 2 which is one of the main results 
of this paper: causality violation and plasma instabilities may arise in the interior of the 
geometry. These are thoroughly discussed in section 5. They should be related to analogous 
unhealthy properties arising in the CFT plasma far from the UV, whose identification and 
full characterization is left as an open problem. They should correspond, however, to very 
different physical situations: while we may expect that causality violation corresponds 
to some illness of the boundary theory, instabilities should be related to either plasma 
instabilities or to issues such as not being expanding around the true vacuum 3 . Anyway, 
interestingly enough, this shows that causality violation is not necessarily a phenomenon 
arising in the UV and it is not generically dual to positivity of energy in the CFT. 

2 This phenomenon already takes place for the cubic Lovelock theory. It leads to a less stringent (thus 
irrelevant) constraint in d = 7, but is crucial to prevent arbitrarily small values of 77/ s in 8 < d < 10 
theories. When d > 10, causality at the boundary is enough to discard this possibility in cubic Lovelock 
gravity [19]. 

3 We thank Rob Myers for clarifying comments on these points. 
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The attentive reader has already realized that higher order Lovelock gravities imply 
higher-dimensional theories. There is a second more subtle point that should be presently 
raised. The analysis performed in this article deals with higher order curvature correc- 
tions whose coupling constants do not need to be neccesarily small. Both conditions are 
unnatural within the string theory realm. It seems however worth inspecting thoroughly 
these theories based on the very fact that they lead to reasonable results. That is, even 
for these presumably non-stringy setups we obtain what seem to be splendid checks of the 
AdS/CFT correspondence. 

Section 6 is devoted to the analysis of the r]/s ratio in higher order Lovelock theories. 
The first part of that section presents analytic and numerical results for the minimum 
ratio in any dimension up to quartic order theory. These are obtained by considering the 
full stability and causality constraints. Our results are relatively smooth functions of the 
spacetime dimensionality. In the second part we consider a limited analysis of arbitrary 
Lovelock theories in any dimension. Picking a specific curve through parameter space, the 
properties of which are carefully discussed in Appendix B, we show that there seems not 
to be a dimension independent viscosity bound. Within our framework, we argue that a 
strongly coupled ideal fluid could be achieved in the strict limit d — > oo, when a corre- 
spondingly infinite number of Lovelock terms are taken into account. This is somehow 
reminiscent of the so-called species problem, though it is not clear to us whether there is 
a more rigorous way in which these two setups are related. Following a discussion and 
conclusions presented in section 7, we have a series of appendices where we derive a holo- 
graphic formula for the parameters entering the 3-point stress energy tensor correlations 
in arbitrary dimensions, and we present the details of the analysis of a particular curve 
through parameter space in an arbitrary Lovelock theory. 



2. Lovelock gravity 

Lovelock gravity is the most general second order gravity theory in higher-dimensional 
spacetimes, and it is free of ghosts when expanding about flat space [17, 20]. The action 
can be written in terms of differential forms 4 as 

K 



^E^K), (2.1, 

fc=0 

where is a set of couplings with length dimensions L 2 ^ k ~ l \ L being some length scale 
related to the cosmological constant, and d is the spacetime dimensionality, if is a positive 
integer restricted to K < [^y^] ■ lZ( k ' is the exterior product of k curvature 2- forms with 
the required number of vierbeins to construct a d-form, 

U {k) = e fl ... fd R hh A • • • A Rhk-ihk A e hh+vfd s (2.2) 



4 We make use of the language of tensorial forms: out of the vierbein, e a , and spin connection, uj a b , 
1-forms, we construct the Riemann curvature, R a b , and torsion, T a , 2-forms: 

R\ = du\ + u ' c A U \ = ii?V dx- A dx v , T a = de a + u\ A e b . 
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where e* 1 '"-" 1 is a short notation for e^ 1 A . . . A e^ k . For convenience, we will also use from 
now on the expression Rj 1 '"^ = Rf 1 ? 2 A • • • A Rh>=-if2k_ The ze roth and first term in (2.1) 
correspond, respectively, to the cosmological term and the Einstein-Hilbert action. These 
are particular cases of the Lovelock theory. It is fairly easy to see that ao = L~ 2 and a\ = 1 
correspond to the usual normalization of these terms. 5 

If we consider this action in the first order formalism, we have two equations of motion, 
for the connection 1-form and for the vierbein. If we vary the action with respect to the 
connection the resulting equation is proportional to the torsion, and so we can safely set 
it to zero as usual, allowing us to compare our results with those coming from the usual 
tensorial formalism based on the metric. The second equation of motion is obtained by 
varying the action with respect to the vierbein. It can be cast into the form 

£a = e afl ... fd : Fffi A • • • A ^( 2 K K f lhK A e/^+r-/,-: = ; (2 . 3) 

where J-^ = R ab — Aj e ab , which makes manifest that, in principle, this theory admits K 
constant curvature vacuum solutions, 

jrab = R ab _ A . g a6 = Q ^ 

The K different cosmological constants are the solutions of the K-th order polynomial 

K K 

T[A] = ^a fc A fc = a K H(A-A 4 )=0, (2.5) 

fc=0 i=l 

each one corresponding to a 1 vacuum'' which we may think of doing perturbations about. 
The theory will have degenerate behavior and possibly symmetry enhancement whenever 
two or more effective cosmological constants coincide. 

2.1 AdS vacua and black hole solutions 

Here we are interested in AdS vacua and black brane solutions. These can be obtained by 
means of an ansatz of the form [21, 22, 23] 6 

ds 2 = -N 2 # f{r)dt 2 + +-^2 fix 2 , x= (X2,...,x d _i) . (2.6) 
A natural frame for the vierbein is then 



N#y/f(r)dt, e l = -±=dr, e a = j dx a , (2.7) 
where a = 2, . . . ,d— 1, and the only non- vanishing components of the spin connection read 

u\ = ^P=e\ (2.8) 



J We set, without any loss of generality, 167r(d — 3)! Gjv = 1 to simplify our expressions (Gn being the 
Newton constant in d spacetime dimensions). In tensorial notation, thus, the cosmological constant has the 
customary negative value 2A = — (d — l)(d — 2)/L 2 . 

6 We can consider different functions in the timelike and radial directions but they are bound to be equal 
by the equations of motion for generic values of the couplings. 
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The Riemann 2-form is given by 

R 01 = _fp. e Ag l i? 0a = _ : f(O e Ae a 

i? la = e 1 Ae a , i? ah = e a A e fe . (2.9) 

2r r 

If we insert these expressions into the equations of motion, we get after some manipulations 



d - 1 /K 



+ (d-l) 



dlogr 

where g = —f/r 2 . This can be easily solved as 



5> fc / =0, (2.10) 



A:=0 



n 9 ]^jz a ^ k = ^{ r ^) d ' 1 , (2.H) 



L 2 

k=0 



where r + is an integration constant that determines the mass of the spacetime. If r + = 0, 
the equation is solved by the constant function g(r) = A, with A a solution to (2.5), and 
the metric reduces to pure (A)dS. The full metric is fixed by choosing N^. = — (L 2 A) _1 , 
so that the speed of light in the boundary theory is set to one. 

For r+ > the equation (2.11) admits a non-trivial black hole solution with a horizon 
at r = r_|_, where g(r) vanishes. Notice that (2.11) leads to K different branches, each 
one associated to each of the cosmological constants (2.5). Just one branch has a horizon 
there. This can be seen since the polynomial root g = has multiplicity one at r = r + . In 
order to have higher multiplicity it would be necessary that a\ vanishes, but a\ = 1 in our 
treatment (it corresponds to the Einstein-Hilbert term). This unique branch of solutions 
with a horizon - from now on the 'relevant branch' - is just a deformation of the usual 
Einstein-Hilbert AdS black hole solution and shares with it a lot of nice features. For 
instance, it has a well defined temperature 

* 4vr Ait L 2 V ' 

as rg'(r) = —(d — 1) ^t^j , and its entropy still satisfies the usual area law [24, 25, 26, 27, 
28, 29] 

r d-2 

S =^p, (2.13) 
where we have divided by the volume of the flat directions. The mass can be obtained as 

This is general for any Lovelock black hole of the type considered here. 

A detailed study of black hole solutions in Lovelock theory is far from trivial (see, for 
instance, [30]). For cubic (or higher order) Lovelock gravities, there are entire regions of 
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the parameter space where there are no black holes at all even if there is a well defined AdS 
vacuum [18, 19]. This is a novel feature that does not take place in Gauss-Bonnet gravity. 
It has to do with the several branches mentioned above. The full solution is not valid 
whenever the relevant branch degenerates with any other one for a given radius. A careful 
analysis leads to the characterization of a non-compact subset of the space of Lovelock 
couplings that do not admit black hole solutions [19]. 

3. Holographic dictionary 

Lovelock theory has a rich structure of AdS vacua. Correspondingly, the AdS / CFT corre- 
spondence tells us that there should be an analogously complex structure of dual CFTs. 
We know very little about these higher-dimensional CFTs. Indeed, at least for the case 
of supersymmetric CFTs, one naively expects to have at most six dimensional non-trivial 
unitary conformal field theories, based on the algebraic construction by Nahm [31]. For 
higher-dimensional CFTs, the general expectation is that they are not lagrangian theories 
and their constituent degrees of freedom are not gauge fields but possibly self-dual p-forms 
(in the case of even spacetime dimensions) [32] . These theories should have a stress tensor 
and conformal symmetry will constrain their 2- and 3-point functions as we shall presently 
describe. 

3.1 Calculation of Ct 

The leading singularity of the 2-point function for a general CFT in any number of dimen- 
sions is fully characterized by the central charge Ct [33] as 

<T Qb (x) 7^(0)) = X a6 , cd (x) , (3.1) 

where 

lab,cd( X ) = \ Moc(x) iftd(x) + Iad(x) i&c(x) - j^— j- T/ab Vcd^j , (3.2) 

and 

J a &(x) = r] ab -2 . (3.3) 

x z 

This structure is completely dictated by conformal symmetry. 

The holographic computation of Ct was performed in [16] for Gauss-Bonnet theory in 
various dimensions. If one considers a general Lovelock theory, the result is very similar as 
we shall see shortly. Since Ct appears in the general 2-point function of the stress tensor, 
it is sufficient to consider a particular set of components. Following [16] we consider the 
correlator {T xy (x) T xy (0)). According to the AdS/CFT dictionary, such a correlator can be 
computed by considering gravitational fluctuations around the AdS vacuum of the theory. 
In the present case, it is sufficient to take a metric fluctuation 0(r, x) = h x v (r,x) about 
empty AdS. Expanding the lagrangian (2.1) to quadratic order in the fluctuation and 
evaluating it on-shell, we find 

1 



.'[A] f dxdtr d {^d r cj>) , (3.4) 
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where Lp 1 2 = 8ttGn is the Planck length and 

K 

T'[A] = Y J ka k A k - 1 . (3.5) 
fc=i 

In the notation of [16] this expression is given by 7 (1 — 2f O0 \ + ...). Indeed, the preceding 
formula exactly matches equation (3.12) in that reference with T'[A] corresponding to the 
specific case of Gauss-Bonnet. Accordingly, our calculation of Ct proceeds in an analogous 
fashion and we find 

_ d T[d] T'[A] 

Ct ~ d=2^ r [*=i] i-^ ' ( 6) 

If there is a dual conformal field theory, or class of theories which are holographically 
described by the Lovelock gravity lagrangian, then the above expression provides the value 
of the central charge for those theories. Unitarity then demands that Ct must be positive, 
and this nicely translates into the gravity side to the condition T'[A] > 0. From the above, 
we see that if this is not so, the kinetic term for the graviton modes acquires a negative sign, 
and hence these modes are no longer unitary and the theory shows the kind of instabilty 
first found by Boulware and Deser (BD) [34] in Gauss-Bonnet gravity. The positivity of 
Y'[A] has been proven to be true for the relevant branch of solutions in any Lovelock 
theory [19]. Thus, this branch is protected from BD instabilities and the corresponding 
dual theory is unitary. 

3.2 Three-point function 

The form of the 3-point function of the stress-tensor in a (d — l)-dimensional conformal 
field theory is highly constrained. In [33, 35], it was shown that it can always be written 
in the form 

(at {1) +bx {2) +cx (3) "l 

{T ab x T cd y T ef z = ^ ' ,w_ ' 3.7 

(|x - y| |y - z| |z - x|) d 1 

where the form of the tensor structures :cd e j will be irrelevant for us here. Energy 
conservation also implies a relation between the central charge Ct appearing in the 2-point 
function, and the parameters A,B,C, namely 

vr^ {d-2){d+l)A-2B-UC 

Since we have already computed Ct in the previous section, we are left with two indepen- 
dent parameters to be calculated. 

A convenient parameterization of the 3-point function of the stress-tensor was intro- 
duced in [11]. One considers a localized insertion of the form f dui ejk e~ tu3t T jfc (x), and 



7 We have the equivalence /oo = —AL 2 
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measures the energy flux at lightlike future infinity along a certain direction n. The final 
answer for the energy flux £ (n) takes the form [15, 16] 



<£(n)> 



n 



d-3 



■ n i n 3 e *ik e ik l_\ + ^ fining] 2 



e* k e ik d-2) \ e* k e ik d{d - 2) 



(3.9) 



with E the total energy of the insertion, and £ld-3 the volume of a unit (d — 3)-sphere. 

The existence of minus signs in the above expression leads to interesting constraints 
on the parameters £2 and £4, by demanding that the measured energy flux should be 
positive for any direction n and polarization e^. In spite of the fact that the positivity 
of the energy flux is not self-evident, there is a field theoretic argument supporting this 
claim [13]. Furthermore, interestingly enough, it was recently argued that this condition is 
equivalent to unitarity of the corresponding CFT [36]. By picking a particular direction and 
considering polarizations with different rotational transformation properties with respect 
to it, one finds three different constraints, 

Tensor : 1 - dh t2 ~ d(^2) u * ' 

Vector: 1 + W^) t2 -W^) u - ^ (3 - 10) 

d-4 d(d-3)-2 
Scalar : 1 + -^r£ 2 + ' U > . 

a — 2 d{d — 2) 

For conformal field theories with a weakly curved gravitational dual, it is possible 
to compute £2 and £4 holographically. This was first done in detail for Gauss-Bonnet 
theory in various dimensions [16] and later applied to quasi-topological gravity in [37]. The 
calculation proceeds by considering the vacuum AdS solution perturbed by a shockwave, 

which corresponds holographically to a T insertion. By adding a transverse metric 

fluctuation, one reads off the interaction vertex from the action, and from that one obtains 
£2 and £4. 

Shockwave backgrounds were also considered by Hofman [13] in the context of five- 
dimensional Gauss-Bonnet theory. There it was found that in the presence of the shockwave 
there is the possibility for causality violation in the dual field theory. This places bounds 
on the parameters of the theory, which precisely match those portrayed in (3.10). This 
shockwave calculation was extended to arbitrary higher-dimensional Gauss-Bonnet gravity 
in [15], and the corresponding constraints were found. In [16], £2 and £4 were computed 
explicitly for these theories and the bounds (3.10) perfectly agree with the constraints 
obtained in [15]. 

We would now like to generalize this story to higher Lovelock theories. We consider, 
along the lines of [15], a helicity two perturbation 4>(u,v,r) in the shockwave background 

N 2 L 2 

ds 2 AdSjSW = * 2 (-dudv + dx { dx i + 2ecj)dx 2 dx 3 + dz 2 ) + f(u) w(x a , z) du 2 , (3.11) 

where z = L 2 jr is the Poincare coordinate and u, v = x° ± x d ~ l are light-cone coordinates. 
This amounts to choosing just one non-vanishing component of the polarization tensor, 
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e 23 7^ 0. Leading contributions (in the high momentum limit) to the equations of motion 
come from the exterior derivative of the perturbation of the spin connection. 8 We get, 



d(5uj 



12 i 



# 



e z 



dl<f> e L A e 2 + I d u d v 4> + f(u) w(x\ z) <9> I e° A e 2 



# 



+ ^2 /(«) ro(x a , *) I e 1 A + (• • • ) e° A e 2 



the ellipsis being used in the second expression since the corresponding term does not 
contribute to the equations of motion. The components with index 3 instead of 2 are the 
only remaining non-vanishing ones, and they are obtained just by changing 9 <\> — > —(p. The 
other ingredient we need is the curvature 2-form of the background metric, that can be 
written as 

R ab = A(e a Ae b + f{u)X ab ) , (3.12) 

where A = — L ^~ N 2 an d X ab is an antisymmetric 2-form accounting for the contribution of 

# 

the shock wave. The relevant component here is 

~2 



X 



la 



(2m + zd z w - z 2 d 2 a w) e° Ae a + (■■■) e° A e z + {■■■) e A 



X l(d-1) = _ 



LN._ 



# 



(3zd z w + z 2 d 2 w) e° A e d ~ x + (•••) e u A e 



a „a 



with a, b 7^ 0, 1, d — 1, and b / a. The relevant equation of motion is given by SS3 A e 3 = 0, 
where 

K 

5£ 3 A e 3 = ka k £3fi-f d -A 6ujhh ) A R h "' hk A e W"/d-i3 



fc=i 



<- 1 )"~ 1 i^E*« Ai - 1 



# fc=i 



(d-sy.idud^ + j^f^zud 2 ^ 



# 



- (fc - l)(d - 5)! -^2" (" 4ro " 2 ^ ro + * + ^tu)) <9 2 < 

We made use of the equation of motion for the shock wave profile 

2(d - 3)w + (d- Q)zd z w - z 2 (d a d a w + d 2 z w) = 0, a = 2, . . . d - 2 . 
This equation admits several solutions of the type 



W = CKQ 



(z 2 + (a; 2 - xlf + • • • + (x n - x%)) 



13 ' 



. (3.13) 



(3.14) 



(3.15) 



8 The 3-point function will be determined by these kinds of terms; actually, terms in d%4> since the vertex 
has to be of the form (f>d^,<f>didjvj, as vo ~ h uu . 



4 dx 2 -> -dx A , dx A dx , 



is a symmetry of the background as well as of the vierbein basis. 
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where n is the number of transverse coordinates. The relevant solution for our discussion 
here is the one given by a = d — 3 and (5 = d — 2. This shockwave profile has been argued 
[11] to be the dual field configuration to £(n) provided x l = 1+ ^_ 2 and f(u) = S(u). 

From (3.13), we shall focus on those terms proportional to d 2 (j). The shockwave inter- 
action term thus gives the following contribution to the equations of motion, 



A T'TAl 

T '[ A ] w + (d- 3 )(d-4) ( 4ro + 2zdzW ~ z2 ^ W + 
There are extra non vanishing components of the form 

5£i A e 2 ~ dzdiwdl4> , 



. (3.16) 



$£ d _i A e 2 ~ -lj (d 3 w + x d ^d 3 d d ^w) d 2 v 



However, these are irrelevant for computing the three-point function we are interested in. 
The tensor channel mixes with other modes but they would only affect other correlators 
irrelevant for our discussion. 

The 3-point function follows from evaluating the effective action for the field <j) on-shell, 
on a particular solution which depends on all coordinates, including X2,x 3 [16]. The cubic 
interaction vertex of (f> with the shockwave appearing in the action will be essentially the 
one in the equation of motion determined above. Up to an overall factor, the cubic vertex 
is then 

5« ~ C T J (l - A T ^ ) (d _ 3 ^_ 4) ) ■ (3.17) 

where 

T2 = z 2 {d\w + dim) - 2zd z w - Aw 
w 

This is nothing but the same T2 appearing in [16]. Indeed, following that paper the relevant 
graviton profile is 

<t>{u = 0, v, x a , z) ~ e~ iEv 5 d - 3 {x a ) 5(z - 1) , (3.19) 
so that we need to impose x a = and z = 1 yielding 

T 2 = 2(d -l)(d- 2) - -L_) , (3.20) 

and we therefore read off 

_ 2(d-l)(d-2)Ar[A] _ 
t2 ~~ (d-3)(d-4) T'[A] ' t4 ~°- (3 - 21) 

This expression reproduces previous results in Gauss-Bonnet gravity [15, 16], and is exactly 
the same as conjectured by de Boer, Kulaxizi and Parnachev in [18]. This completely 
proves the equivalence between the positivity of energy constraints and those coming from 
causality in the boundary theory (see next Section). Using these results together with the 
expression for Ct, we find expressions for the usual 3-point function parameters A,B,C in 
Appendix A. 
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4. Shear viscosity and finite temperature instabilities 



It has been shown in [38, 39] that the shear viscosity to entropy ratio of a generic Lovelock 
theory is given by 

where A = (I2/L 2 , and we have defined the critical value A c where the ratio vanishes 

It is interesting to stress that the ratio only depends on the quadratic Lovelock coupling, 
02- Generically one is not free to take an arbitrary A. There are causality constraints which 
restrict the range of values where the theory is well-defined [4, 7] . These restrictions match 
the restrictions on t2 and £4, and these necessarily involve all the parameters in the theory. 
The maximum value A( ) that may be achieved depends on the theory in question. In this 
way it has been found [18, 19] that causality constraints alone cannot prevent one from 
reaching (and even overpassing) r]/s = 0. In this section we show that as one approaches 
small values of rj/s, the black hole becomes unstable and the linear approximation breaks 
down. This places an effective set of constraints in the parameters of the theory, saving 
the ratio above from ever becoming too small. But first, let us briefly recall how to derive 
the constraints by demanding causality. 

4.1 Graviton potentials and causality 

We start off with the black hole background (2.6), and add metric fluctuations h a b of 
frequency uj and momentum q in a fixed direction. The fluctuations split into three channels 
according to their polarization relative to the momentum, namely the tensor, shear and 
sound channels [40] (or equivalently helicity/spin two, one or zero). The equations of 
motion for these dynamical degrees of freedom, 4>i(r), the subindex i = 0,1,2 indicating 
the corresponding helicity, can be recast as Schrodinger type equations [41]. In the large 
momentum limit they reduce to a very simple form 10 

-??&* i + U i (y)* i = a?*i, n=-^0, (4.3) 

y q 

where a = ui 2 /q 2 , y is a dimensionless tortoise coordinate defined as dy/dr = \J — A//(r), 
and = Bi(y) 4>i(y), where Bi(y) = Bi(g(y)) are functions of the metric whose specific 

expression can be found in [41]. In order for the linear analysis to be valid we need a 
non- vanishing function Bi(y); we shall comment on this later on. The effective potentials 
Ui, can be determined as the speed of large momentum gravitons in constant y slices, 

V < ° ' (4-4) 

y = , 




10 In the notation of [41], we must identify ji = 7 = q 2 L 2 , and our potentials are related to theirs, 
Ui — > — Vi/7 A, as 7 — > 00. 
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y = being the boundary of the spacetime. 

For a generic Lovelock theory, in terms of the original radial variable, one finds for the 
tensor, shear and sound channels, respectively [19]: 



(d-3) Cf } [g,r] 
r 2 M = 9/ A ( 2C^[g,r] _ cf\g,r] \ 

where we have defined functionals, C^fg, r], involving up to feth-order derivatives of g: 

Cf[g,r) = r[g], (4.6) 

dpM = { r t + {d ~ 2 ' ) ) r[9] ' (47) 

cf M = ( r J; + (d-3)) (r^ + (i-4))rM . (4.8) 

Notice now that we can make a change of variables (r,g) — > (x = logL 2 T,i ? = logL 2 T') 
so that rd r = — {d — l)d x (as r and T are related by (2.11)) yielding for the potentials 



4(x) 



(d - 1) g (g - ng) (g - ng) + 

(d-4)A ^_ F / (x) ) 



( d - 1) g [g ~ ng J (£f - gOg ) - F"(X) 

(d-2)A 



In order to avert causality violation, we must demand these effective potentials to be always 
smaller than one [4, 7]. In particular, at the boundary we have c 2 = 1 and there we must 
demand 11 d x cf < 0, as x — > — oo. Using d x g = T[g\/T'[g\, we get the following constraints: 

Vector : T'[A] - AT" [A] > . (4.10) 

(a — 3) 

Scalar : T'[A] - jj - ^ AT" [A] > . 
(a — 3) 



1 Recall that the boundary corresponds to Xb = — oo and the horizon to Xh = 0. 
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These can be rewritten in terms of the dual CFT parameters, using the expressions for 
£2,^4 in (3.21). The constraints become: 



Tensor : 1 — ti > . 

d 2 

Vector: 1 + t 2 > . (4.11) 

d — 4 

Scalar : 1 + ti > . 

d 2 

These precisely match the constraints (3.10) coming from positivity of energy in the dual 
conformal field theory (with t± = 0). 

4.2 Plasma instabilities 

It has been found in [42, 43, 16] that, for certain values of the Gauss-Bonnet coupling, 
some effective potentials might develop negative values close to the horizon. This has been 
previously argued to indicate an instability of the plasma [44]. In eq.(4.3) we see that the 
role of h is played by 1/q. By taking sufficiently large spatial momentum (small ft), we can 
make an infinitesimally small (negative energy) well to support a negative energy state in 
the effective Schrodinger problem. Going back to the original fields, this translates into an 
exponentially growing and therefore unstable mode [44]. 

We now show that the same phenomenon occurs in a general Lovelock theory 12 . A 
negative potential well develops whenever the slope of the effective potential at the horizon 
is negative, and so we must require d x cf < there. At the horizon we have g = 0, and it is 
straightforward to make use of T(g) = L~ 2 e x to find the derivatives of g(x) at the horizon, 
e.g., d x g\ g= o = (T[g]/T'[g]) \ g= o = clq/cli = L~ 2 . In this way, we can relate the values of 
the derivatives of F{x) at x = (the horizon) with the coefficients of the polynomial T(g). 
In particular, 

F'(0) = 2A, F"(0) = 2(/i + A(l - 4A)) , (4.12) 

where, again, A = 02 /L 2 , and we have defined fi = 3as/L 4 . Using these results, we expand 
the graviton effective potentials close to the horizon to get 



i-l x (a-^(fef-2A) + 2(, + A (1 -4A)) 



\d-l 



2A 



2 Some work in this direction has been done recently in [45, 46] 
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Stability requires that the potentials are positive. This gives rise to the following constraints 
for the Lovelock couplings: 



(d-3)-2A(d-l) >0 , (4.14) 
(d - 3)(d - 4) - 2A(d - l)(d- 6) - 4A 2 (d - l) 2 + 2fi(d - l) 2 > , (4.15) 
(d - 2){d- 3) - 6A(d - l)(d- 2) + 12A 2 (d - l) 2 - 2^(d - l) 2 > . (4.16) 

This is one of our main results. These inequalities represent the constraints from the 
shear, tensor and sound channels respectively. They are required to hold, otherwise the 
black brane solution is unstable. For /j, = 0, these results match those for Gauss-Bonnet 
gravity derived in [16]. It is important to notice that, in spite of considering a completely 
general Lovelock theory with as many terms as we wish, these stability constraints involve 
just the lowest two Lovelock couplings, A and /j,. 

The first constraint immediately rules out the region of the parameter space where the 
shear viscosity to entropy ratio becomes negative, since A < A c . It is interesting to point 
out that the shear channel stability condition is also needed to ensure the validity of the 
linear analysis [45], as long as all the Bi{r) functions are seen to be proportional to some 
power of the shear potential. Then, we can also relate the would be negative values of 
the shear viscosity - that would lead to a manifest instability given its role as a damping 
coefficient in the sound channel (a negative value then amplifies the sound mode) -, to the 
accompanying break down of the linear analysis (at the horizon). In addition to this, the 
remaining constraints define a new allowed stability wedge in the parameter space: 



(d - 3)(d - 4) - 2 A(d- l)(rf - 6) - 4A 2 (d - l) 2 
2{d- If 



with apex at 



(d - 2){d - 3) - 6A(d - l)(d - 2) + 12A 2 (d - l) 2 

- 2(d^Tp ' (418) 



d-3 (d-3)(d-5) ,„ . 



or, equivalently, on the intersection of the rj/s = line with \i = A (4 A — 1). Then we 
can approach, from the inside of this stable region, a unique point (in the cubic case; a 
(hyper-)line in higher order cases) where r//s = 0, and this is exactly at the apex of such 
region, where the stability constraints coming from each helicity meet. 

For d = 7, the apex coincides with the point of maximal symmetry (A, fj,) = (1/3, 1/9), 
where the polynomial T[g] has a single maximally degenerate root. This is also the Chern- 
Simons point of the d = 7 Lovelock theory. 13 We plot the different regions in the space 
of parameters in figure 1. The minimum value for r\js in this case is attained when the 



13 The Chern-Simons point is the maximally symmetric point of a K-X)ci order Lovelock theory in d — 2K+1 
dimensions, where it actually becomes a Chern-Simons theory for the AdSd group (see, for example, [47]). 
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Figure 1: The allowed region by causality and stability for cubic Lovelock theory in d = 7. The 
black points are the maximally degenerated point (1/3,1/9) (also the Chern-Simons point in this 
case) and the intersection of the helicity zero stability constraint and the helicity two causality 
constraint ~ (0.20,0.017), that gives the lowest value for rj/s. 

upper (helicity zero) stability curve intersects the lower (helicity two) causality curve, which 
happens at 

A ~ 0.202042 , n ~ 0.0175986 , (4.20) 

giving a minimum value 

- ~ 0.393874 x ~ < 0.4375 x ~ , (4.21) 
s 4tt Air 

where the latter one is the correction to the KSS bound coming from GB alone (this 
corresponds to the value of A at which the lower blue curve intersects the axis, Agb = 3/16). 
This seems to be the end of the story in d = 7, at least in the context of Lovelock gravities. 
The existence of plasma instabilities (they originate in the behavior of the gravitational 
potential close to the horizon) sets a definite (positive) lower bound on cubic Lovelock 
theory. However, as we will see, the situation becomes more involved in higher dimensions, 
where the restrictions discussed so far are not enough to prevent rj/s from becoming zero. 

5. Bulk causality and stability 

5.1 The cubic theory in higher dimensions 

Our results imply that the rj/s ratio can never be too small for cubic Lovelock theory in 
d = 7 by a combination of two constraints - preservation of causality, holographically dual 
to positivity of energy; and stability of the black brane solution which can be broadly 
identified with the stability of the thermal plasma. Of course the obvious question is 
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whether similar results are valid in higher order theories. As it turns out, these two 
constraints are indeed enough, but not in the simple fashion we have described so far - 
we must require causality and stability everywhere in the bulk, and not just from a simple 
near boundary or near horizon analysis. 




0.15 0.20 0.25 0.30 0.35 0.40 0.28 0.30 0.32 0.34 0.36 



Figure 2: Causality and stability regions in d — 9 cubic Lovelock theory (and zoom). The thick 
lines correspond to causality at the boundary and stability at the horizon, respectively in blue and 
red, while the dashed lines correspond to causality (in the tensor channel) and stability (in the 
sound channel) in the full geometry. 




_^ , , , , , , , r — , r , , , , , , , , , , , , ^_ A _j , r , , , , , , , , , , , ^_ A 
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Figure 3: Causality and stability regions in d = 11 cubic Lovelock theory (and zoom). The thick 
lines correspond to causality at the boundary and stability at the horizon, respectively in blue and 
red, while the dashed lines correspond to causality (in the tensor channel) and stability (in the 
sound channel) in the full geometry. 

This can be seen already in the simplest case of cubic Lovelock theory in higher di- 
mensions. In figures 2 and 3 we show the relevant part of the causal and stable regions 
for d = 9 and d = 11, as determined by the near horizon and near boundary expansions of 
the effective potentials. There are also a few dashed lines plotted whose meaning we will 
explain in a moment. It is apparent that by going to higher dimensions, a disconnected 
region close to the apex of the stability region is now causal as well. This means that such 
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a region satisfies all our previous constraints, and yet it will have an arbitrarily low shear 
viscosity. This is because, as we have seen earlier, the apex of the stability region actually 
lives on the surface r/ = determined by the stability constraint in the shear channel. This 
point lies inside the causality allowed region for d < 10 [19]. 

While in principle there is nothing wrong with having a small ratio, this is not the end 
of the story. As we have mentioned before, by considering the full effective potentials we 
can see that they develop causality problems and/or instabilities, but now in the interior 
of the black hole geometry (see figures 4). In the left figure, even though the potential 




-0.2 L -1.4 -1.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 



Figure 4: These figures illustrate the appearance of bulk instabilities (on the left) and bulk 
causality violation (on the right) for the case of cubic Lovelock theory in ten dimensions. They 
correspond, respectively, to the sound potential for A = 0.33 and fi = 0.108 and to the tensor 
potential for A = 0.35 and fi = 0.25. The plot of the potential ends at the boundary where the 
function g takes the value of the cosmological constant, respectively A = —2.5 and A = —1.512. 

initially rises close to the horizon [g = 0), it dips into negative values. In the same way, 
in the figure on the right, the potential goes below one close to the boundary and then 
makes a hump above this value, leading to causality violation. Such features cannot be 
fully seen in a perturbative analysis, although their presence is easily guessed at. Consider 
for instance moving along the curve where the causality bound is saturated in the tensor 
channel. This means that the effective potential is of the form 

V eff = l + 0xe x + v 4 e 2x + ... (5.1) 

It is clear then that along this curve we will see causality violation if ever becomes 
positive. Beyond the point in parameter space where this occurs, there will be a competition 
between the e x and e 2x terms in the potential expansion which will determine the shape of 
the causality curve from that point onwards. In general it is clear that to determine this 
curve one must consider the full effective potential, since all terms in the near-boundary 
expansion become of the same order. Analogous statements hold for the instability analysis 
near the horizon. These curves are determined numerically and shown as dashed lines in 
figures 2 and 3. We can see, in particular, that the sound mode instability curve neatly 
cuts off the region where rj/s might become too small, thereby imposing an effective lower 
bound on this ratio. 
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Going to higher order Lovelock theories, one may wonder if this general reasoning will 
hold. Could the addition of further couplings allow us to escape even the full causality and 
stability bounds? This does not seem to be so. Evidence for this conclusion comes from 
an analysis of quartic Lovelock theory, to which we now turn. 

5.2 Quartic Lovelock theory 

Quartic theory corresponds to the introduction of an extra non-zero parameter, 04 = ^L 6 . 
The theory exists in d > 9, and examining it requires solving for quartic polynomials. The 
analysis of the parameter space of this theory is made in a somewhat similar manner to the 
cubic case, and we will not go into great detail here. There are resemblances and differences 
to the cubic theory. For instance, as in the lower order cases the boundary of the excluded 
region is just the surface where the cosmological constant associated with the relevant 
branch becomes degenerate. However, this surface ends at one of the triply degenerated 
lines (three cosmological constants equal) and we need another surface to actually close the 
excluded region. This also happens in the cubic case where the excluded region is closed 
by /i = A 2 , which is the line along which all three solutions become degenerate at some 
point inside the geometry. In the quartic case the situation is analogous and the excluded 




Figure 5: Surfaces delimiting the (boundary) causal region for the helicity zero and two modes in 
d = 9 dimensions; respectively t 2 = —7/5 (blue) and t 2 = 7 (green). We omit the excluded region 
boundary surface for the sake of clarity. 

region is closed by the surface swept by the line where three solutions degenerate at some 
value of the radius. This surface cuts off the other surfaces bounding the causal regions for 
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the different helicities: green for tensor channel and dark blue for the sound. As always, 
the causal region is contained between the tensor and sound surfaces (figure 5). Besides, 
we also realize that the values of A are not bounded (neither from below nor from above) 
just due to causality - by increasing v it is possible to go to arbitrarily high A. 

This situation drastically changes once we consider stability and causality constraints 
(see figure 6). When horizon stability is imposed, the allowed region of parameters is 




Figure 6: Stability close to the horizon (in purple) and (boundary) causality surfaces for d = 9. 
It can be observed that there is a channel connecting the two disconnected regions found in cubic 
Lovelock theory (see figure 2). 

further constrained by the purple surfaces; notice, in particular, that it includes a channel 
connecting the Einstein-Hilbert point with the rj/s = line (the apex of the stability region, 
in purple). We are therefore in a situation analogous to the one in the cubic theory, and 
we expect the full stability and causality constraints to rescue us. This is indeed the case. 
In figure 7 we show the causality and stability constraints taken at different v slicings. As 
v increases, the causality curves move up in the plot until at some point they block the 
path to the critical point where rj/s = 0. For lower values of u, instead, we cannot reach 
this point due to the full sound stability constraint. 

We have actually done this analysis for various v and have reached the same conclu- 
sions. These results add further evidence that for every dimension, stability and causality 
place an effective lower bound on the value of the shear viscosity to entropy ratio. 
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Figure 7: Different slices for the quartic case in 9d. Top left corresponds to v = —0.02, top right 
v = and below v — 1/64. The black thin dashed line is an imaginary line connecting the different 
end points of the causality constraints. The thick black line is the boundary of the excluded region. 
The dashed red and blue lines are the full sound stability and tensor causality curves respectively. 

5.3 Expansions at r]/s = 

We would now like to present a general argument that shows one should expect instabilities 
to appear generically as the r//s = point is approached. The argument is related to the 
fact that in the section 4, there is something which might invalidate the analysis. This is 
when the perturbative expansion on the radius, close to the horizon, breaks down. Notice 
that the tensor and sound mode effective potentials, c| and Cg in (4.9) have (d — 3)/(d — 
1) — F'(x) in the denominator. At the horizon we have, 



This means that the near-horizon expansion breaks down whenever A ~ A c , precisely where 
rj/s ~ 0. Close enough to this point the linear analysis breaks down, and we must treat this 
case separately. This is true for any Lovelock theory in any dimension, and it is therefore 
of interest to see if there are any general statements one can make. 

Let us then investigate the behavior of the effective potentials given in (4.5) at the 
horizon. First notice that the equation (2.11) fixes 



d-3 
d- 1 



F'{x) ~ 2(A C - A) . 





(5.3) 



and we can expand 



(d-3 
\d-l 




d-3 
d-l 




71=1 



X 



,n 



(5.4) 
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The expansion is such that each term in the above is controlled by a different Lovelock 
parameter. For instance -F'(O) = 2A, and so the first term vanishes whenever A = A c . When 
the following derivatives F^(0) vanish, Vj < n, we can actually write 

p(n) (Q) = T (n+1) [Q] _ "Q d ~^_~ l , (5.5) 
i=l 

and so we can set the following term to zero by choosing 

X(»)[0] 1 d-2i-l . . 

sr-sn-sri-- < 5 - 6 > 

i=l 

In particular, 

_ l (d-3)(d-5) _ l (d-3)(d-5)(d-7) 

^"2 (d-1) 2 ' ^ c_ 6 (d-l)3 ' ■■" 1 J 

For a given (odd) dimensional Lovelock theory, if we set all K = couplings to these 
particular values, we will be at the so-called AdS Chern-Simons point [47]. Now, notice 
that (5.4) is actually the denominator of the helicity two and helicity zero potentials (and 
proportional to that of the helicity one mode). The leading correction to those potentials 
close to the horizon, when we set F'(0) = 2A C and F^' <n )(0) = 0, is then 

o , ^ d 1 Tt 1 - , , . 

c ^ x ) ~ ~jzr A + 0{x) ' (5 ' 8) 

»'-^3 (,-l)lpA ' , " + °W»' (5 ' 9 » 
c « {x)a ^lJX + 0{x) - <5 ' 10) 

The effective potentials have therefore discontinuous limits at the horizon as one approaches 
A = A c . Since A < 0, the sound channel potential tends to a negative constant at the horizon 
as the critical value A c is approached 14 . We conclude that for any given Lovelock theory, 
it is impossible to get arbitrarily close to rj/s = without running into an instability. This 
means that, at least within Lovelock theories of gravity, it seems to be impossible to obtain 
arbitrarily small values for the rj/s ratio at any given dimensionality. 

These arguments say nothing about the value of the minimum itself. It seems clear 
that for any Lovelock theory there will be a minimum, but it is also reasonable to expect 
that for higher dimensionalities, the existence of extra free parameters would allow a lower 
value to be reached. For instance, looking naively into the horizon expansion of the effective 
potential 

V e ff = ax + fix 2 + . . . (5.11) 



14 A caveat to this conclusion is the case where one takes the full function F'(x) — (d — 3)/(d — 1) to zero, 
what actually corresponds to the Chern-Simons point. However, it is easy to see from their definition that 
in this case one of the c| and Cq potentials becomes negative. 
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it seems that one could move towards higher values of A staying inside the stability region by 
simply picking parameters such that a, (3, etc. are always negative. In practice, however, 
we necessarily run into difficulties with this reasoning because of the breakdown of the 
perturbative expansion as we approach the critical value A c . 15 In general, the precise 
interplay between these two issues must be determined by numerics, on a case by case 
basis. 

6. The r]/s ratio in higher order Lovelock theories 
6.1 The minimum ratio in quartic theory 

In Gauss-Bonnet gravity and cubic Lovelock theory, the parameter space is sufficiently 
small that it can be explored easily. In particular, one can find the full numerical causality 
and instability curves which determine the allowed regions in parameter space. As one 
goes to higher orders this analysis becomes more and more baroque, as the hypersurfaces 
in parameter space delimiting the causal and stable regions are more complicated and can- 
not be visualized in general. However, as was clear from section 5.2 we have succeeded 
in studying in detail quartic Lovelock theory in various dimensions. For each theory and 
dimensionality we have found the minimum value for the shear viscosity to entropy ratio, 
and plotted the results in figure 8. The first thing we can notice is the remarkable smooth- 




Figure 8: The dots correspond to the numerical (some analytic) bounds for Gauss-Bonnet (red), 
cubic (green) and quartic (blue) Lovelock theories, and the colored dashed lines correspond to 
the estimated asymptotic value on each case. The dashed black line corresponds to the curve 
Airrj/s = 2/(d — 2) that nicely fits (see the zoomed figure at the right) the first seven absolute 
bounds (including the unit value for d = 4) and has a tantalizing behavior, 4ttt]/s ~ 2/d ~ 1/-Kmax, 
for d — > oo, where K m&yi is the integer part of (d — l)/2. 

ness either for each separated value of K as well as extending from one value K to K + 1. 
There is nothing preventing the actual values to jump sharply given the discreteness in 
the number of couplings involved in these theories. Interestingly, we find that there is a 
minimum as a function of dimension for any given order. The critical dimension for these 

15 In other words, even though the first term in the horizon expansion might iead to a positive effective 
potentiai, the higher order terms might reverse this tendency, and they are becoming more and more 
important as A — > \ c . 
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K 


dmin 


(4irr)/s) min 


2 


9 


219/529 ps 0.414 


3 


14 


0.222 


4 


23 


0.176 



Table 1: Critical dimension (d m i n ) for which the minimum value of the viscosity to entropy density 
ratio is attained for each order in Lovelock (K) theory, as well as the corresponding actual minimum 
value of the quantity (4arr)/s) m i n . 



minima and the corresponding values of the minimum rj/s for each theory are quoted in 
table 1. Overall, the minimum rj/s ratio seems to be decreasing rapidly as one increases 
the number of couplings. 

It is interesting to comment on a very simple function that nicely fits the minimum 
value of rj/s for d < 11. It reads 

s 4vr d - 2 v 1 

It is the simplest curve that smoothly interpolates these points. It also has a nice asymptotic 
behavior, as we will discuss shortly. However, it is important to stress that this is not a 
(dimension dependent) bound for rj/s. Indeed, this is already clear in the most relevant 
case given by d = 5, where 2/3 is actually (slightly) greater than 16/25 (see the right 
zoomed figure 8). The expression in (6.1), though, approximately captures the dependence 
of a dimensional dependent novel bound for rj/s arising in Lovelock theories. 

6.2 No dimension-independent rj/s bound 

In higher order Lovelock theories, the parameter space becomes too large and we must 
resort to a more limited analysis. We have chosen to concentrate on a particular line 
through parameter space, parameterized by the single parameter A: 

T[ 5 ] = l+ 5 + A 5 2 + ... = (l-£) (l-£) 1 . (6.2) 

All other parameters can be determined from A (A is fixed by the overall normalization). 
The curve has the property that it starts at the Einstein-Hilbert point and ends at the 
maximally degenerate point (MDP) of the theory, at A = —K where the polynomial has a 
single degenerate root. A detailed analysis of the stability and causality properties of this 
curve is performed in appendix (B). We have found that for generic values of K and d 
there is an interval of A values where there is instability at the horizon 16 . This is the range 
of values where our curve goes outside the stability wedge. Beyond this interval, our curve 
returns to the stability wedge, but nevertheless there are still instabilities - no longer at 
the horizon but in the bulk. These instabilities persist all the way to the end of the curve, 
where it reaches the maximally degenerate point. 



16 Except for K < 7, where it disappears for high enough d. Whenever the unstable interval is not present 
the whole curve is stable. 
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From these results it is clear that generically, the maximum stable and causal point 
one may reach along this curve is the first crossing point with the stability wedge. This 
crossing point can be found analytically, though the results are not particularly enlighten- 
ing. However, there is a considerable simplification at very large dimensionality. To first 
approximation the crossing point is then determined by solving the algebraic equation, 



K 2 (3 + 2A) + 2KA(4 + 3A) + A 2 (6 + 6A + A 2 ) = 



(6.3) 



Solving this equation, and using (B.3), we can find the value of A at the crossing point. 
However, this value actually decreases as K increases. This leads to a minimum rj/s ratio 
which increases as we include more couplings. Notice there is no contradiction here: recall 
we are moving along a particular curve in parameter space, which is simply not the optimal 
one in terms of finding the minimum possible value for this ratio. 

Although the curve we have chosen is not optimal in this sense, it is still useful as we 
can use it as a base point for exploring nearby regions of parameter space. In particular, 
consider moving along the curve towards the maximally degenerate point. For high enough 
d, it is shown in Appendix B that close to this point there is only an instability in the 
sound channel, i.e. all other channels are causal and stable 17 . We now perturb our curve 
by considering a small perturbation 5Xg 2 or 5 fig 3 to (6.2). The resulting potentials can 
then be easily found numerically. 




Figure 9: On the left, the helicity zero potential for K = 8, d = 100 and A = — 8 (dashed), 
—7.9, —7.5, —7, —6.5, —6, —5 along our curve. On the right the same curve slightly perturbed by a 
5Xg 2 term in the polynomial. The precise values of i5A are 5X = 0, 1.8 • 10~ 18 , 6 • 10~ 13 , 1.9 • 10~ 10 , 6 • 
10- 9 ,8-10- 8 ,3-10" 6 . 



In figure 9 we show the spin-0 effective potential for K = 8, for various A before and 
after adding a small fluctuation. The dashed curve represents the limiting curve for K = 8 
along our particular trajectory in parameter space (A = —8). On the left-hand figures, we 
see how although at the MDP itself the potential does not present any instabilities, it is 
unstable all the way up to that point along our curve, i.e. there are negative potential wells 
in the bulk. On the right-hand side figure we see how adding just a small perturbation is 
capable of lifting the instability. This sensitiveness of the potential to fluctuations is due to 



17 The other relevant constraint is actually the tensor causality one but the causality violating region can 
be made as small as one wishes just by increasing the dimensionality. 
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the fact that the maximally degenerate point is special. At this point we have T^ n )(A) = 0, 
with n = 0, . . . , K — 1, which means that the perturbative boundary expansion of the 
potentials breaks down. Accordingly there is a discontinuous limit at that point. Indeed, 
taking the A — > —K limit along our original curve we find that the effective potentials 
become (assuming d ^ 2K + 1): 

d-3K-l 
° 2 = -tf»(d-4) 9{X) ' 

2 d-2K - 1 . . . JX 

c? = , (6-4) 

Since g(x) € [A, 0] the above are actually always positive for sufficiently large dimensional- 
ity. In particular, although the cq potential always has a negative region for A larger than 
but close to —K, it is actually tending discontinuously to a perfectly reasonable potential. 
This is clear from figure 9. 

Our results then indicate that it is reasonable to expect that sufficiently close to the 
MDP there should be a trajectory in parameter space which is both stable and causal. We 
have shown this is certainly true up to K = 8. If this result holds for higher K, then it 
is always possible to reach sufficiently close to the MDP and have well behaved effective 
potentials. As such, the theory seems to be well defined at that point. 

At the MDP we have 

K-l 

Amdp — 2^ (6-5) 
and therefore a value of the viscosity to entropy ratio 

"q\ 1 1 



(6.6) 

S I MDP 47T K 

By taking high enough K, this can be made arbitrarily small. Besides, since the maximum 
K, -ftT max , is proportional to d/2 in the large d limit, we obtain the asymptotic behavior 
already suggested by the formula (6.1) (see also figure 8). 

To summarize, we have given evidence that it is possible to find a point in parameter 
space, sufficiently close to the MDP of a particular Lovelock theory, such that no stability 
or causality issues occur for high enough d. In particular we have checked that it is possible 
to do this up to K = 8. We conjecture that this can always be achieved for any K, and we 
therefore come to the conclusion that the viscosity to entropy ratio can be made arbitrarily 
small by considering a Lovelock theory of high enough order. In other words, there is no 
bound for the n/s ratio which is independent of the dimensionality, at least in the class of 
Lovelock theories of gravity. This goes in line with our expectations that adding more and 
more couplings, or free parameters in the Lagrangian, it should be possible to reach lower 
and lower values for the shear viscosity to entropy ratio. 
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7. Discussion 



In this paper we have considered various aspects of arbitrary Lovelock theories of gravity, 
both at zero and finite temperature. For these theories, we have explicitly computed the 
parameters t<i and £4 appearing in the M aldacena-Hofman parameterization of the 3-point 
function of the stress tensor [11], and we have checked that they satisfy the positivity of 
energy bounds conjectured in that paper. These results constitute a proof of previous 
conjectures in the literature [18, 19]. 

It is intriguing that classical Lovelock gravity in AdS space can capture highly non- 
trivial aspects of conformal field theory physics. While Einstein gravity is known to describe 
the universal spin-2 sector of a large class of strongly coupled CFTs, no dual to Lovelock 
gravity (beyond GB) has been found, nor have the Lovelock terms with finite coefficients 
appeared in string theory: While higher powers of curvature generically appear in string 
theory as a' corrections to the supergravity action, these corrections are necessarily pertur- 
bative. In Lovelock theory we have dimensionful parameters appearing in the lagrangian, 
which do not appear to be protected from quantum corrections by any symmetry. That 
from the classical theory we can obtain specific numerical predictions on which values of 
the parameters should or not be allowed by causality and that these restrictions should 
agree with analogous restrictions in a hypothetical dual CFTs is highly suggestive. We 
obtain what seem to be splendid checks of the AdS/CFT correspondence, the results being 
smooth functions of the spacetime dimensionality. 

One possible hint is that Lovelock theories all satisfy £4 = 0, which is required for 
supersymmetric CFTs [48]. For free field theories, £4 = is obtained by picking the 
appropriate supersymmetric matter content. In Lovelock gravity, the £2 and £4 parameters 
can be obtained in the boundary expansions of the effective potentials (4.5). The £4 = 
result should imply some relationship between these. Indeed, it is easy to see that from 
(4.5) we have 

(d - 4) c\(r) + {d-2) c 2 (r) - 2(d - 3) c\{r) = , (7.1) 

which is actually valid at any radius r and not just at the boundary, where it indeed 
implies £4 = 0. Beyond the boundary expansion, the relationship written above should 
impose constraints on higher n-point functions of the theory. It is tempting to conjecture 
that these are the holographic duals of some sort of supersymmetric constraint. 

An important part of our analysis was the study of these effective potentials, searching 
for causality violation or instabilities in the large momentum limit which might rule out 
regions of the parameter space. Notice, however, that while the former condition is a 
bona-fide constraint on the theory, since it would lead to causality, the second seems to 
hold only for the validity of the black brane solution, and is not a fundamental restriction. 
We have found that causality violation can occur deep inside the bulk, and therefore this 
effect cannot be captured by a boundary expansion. This means that, contrary to previous 
expectations, causality violation is not necessarily a UV feature of the CFTs related to 
positivity of energy restrictions. Holographically this means that there are interesting 
constraints on the parameters of the CFT which cannot be seen in perturbation theory. 
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The field-theoretic origin of such constraints remains a mystery. 18 On the gravity side, 
these restrictions were essential to prevent an arbitrarily small rj/s ratio for a fixed theory, 
as seen for the concrete examples of cubic and quartic Lovelock theories. 

For higher order theories, we have resorted to the study of a particular curve through 
parameter space. We chose this curve on the basis of simplicity and the capability of 
analytic treatment. The upshot of our analysis is that it seems likely that there is some 
curve leading arbitrarily close to the maximally degenerate point of the theory. The reader 
might be suspicious about this - after all, at this point the central charge of the dual 
conformal field theory is vanishing, and so is the kinetic term of the graviton fluctuations 
around the AdS black hole. That is why we emphasize that one is really considering a curve 
leading to this point, but causality constraints keep us at a safe distance where everything 
is well defined. 

In all cases we considered in detail, the minimum value for the ratio was obtained well 
away from the maximally degenerate point. In any case, we hope our work stimulates study 
of the properties of the theory around this particular point. In the case d = 2K + 1, these 
theories display symmetry enhancement becoming gauge theories of the Chern-Simons 
group: they have local AdS symmetry whereas all other theories with d > 2K + 1 have 
local Lorenz invariance [47]. Contrary to what is suggested by the fact that Ct vanishes 
in the dual conformal field theory, there are some hints pointing towards the existence 
of interesting physics on these theories. For instance, the thermodynamics of their black 
holes displays a qualitative difference to that of generic AdS Lovelock black holes: the 
temperature grows linearly with the horizon radius, and the specific heat is a continuous, 
monotonically increasing and positive function of r+ [49]. Thus, Chern-Simons black holes 
can reach thermal equilibrium with a heat bath at any temperature and they are stable 
under thermal fluctuations. There is a mass gap between the massless black hole and AdS 
spacetime, as it happens in the case of d = 3 [50]. 

Our results indicate that, in the case of very large d, it should be possible to reach a 
value of the rj/s ratio of at least 

2 = JLi (7.2) 

in the Kih. order theory. Intuitively, this is simply telling us that the more parameters we 
have, the lower the ratio can become. This is reminiscent of the species counter-argument 
to a minimum bound on rj/s [51, 52]. Of course, in Lovelock gravity one needs to increase 
the number of dimensions in order to have more parameters. Then it would seem that for 
all practical purposes there is a bound on rj/s for any finite dimensionality. This is certainly 
a possibility, and it is definitely true for Lovelock theory - a combination of stability and 
causality insures it. However, even for fixed dimensionality it seems like it might be possible 
to have interesting higher curvature corrections with finite coefficients [53, 54, 55]. These 
"quasi-topological gravities" have exact black hole solutions of a similar nature to those of 
Lovelock theory, including their thermodynamic features. However, hydrodynamic as well 



18 It would be of definite interest to study whether these constraints are related or not to unitarity 
restrictions that were very recently shown to arise in the computation of 2-point functions of the stress 
energy tensor at finite temperature and large energy and momenta [36]. 
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as causality and stability properties, are quite different [37]. While there is much work 
to be done understanding these theories, it seems plausible that a sufficiently high order 
quasi-topological theory will have a very small rj/s ratio. 

The stability wedge in all Lovelock theories has a non-trivial shape in the A-/i plane 
but is otherwise a cylinder in the remaining couplings. This means, in particular, that 
there is a lower bound for the cubic Lovelock parameter ji. It results from instabilities of 
the black brane solution and, as such, should be related to properties of the CFT plasma 
rather than of the zero temperature CFT. The precise meaning of this bound in the field 
theory side and its consequences in the computation of physical quantities such as transport 
coefficients or thermodynamic variables, as well as the absence of a higher bound for this 
parameter, are interesting avenues for further research. 

As discussed earlier in the case of GB [15] and cubic Lovelock theory [19], it is tempting 
to say something about the existence of a lower bound on A (see figure 10), that implies an 
upper bound for rj/s. It is due to the causality constraint in the former and the stability 
constraint on the latter. This seems to be in line with the expectation that the shear 
viscosity of a strongly coupled system cannot be too large. On this respect, it is interesting 
to mention that this bound disappears in Lovelock theories of quartic or higher order - i.e., 
for d > 9 (see figure 10) -, since no expected problems as causality violation or instabilities 
of the sort discussed in this paper arise in that direction. This might be hinting towards 
the existence of pathologies of a different nature that we have been unable to characterize. 



Figure 10: The dots correspond to the numerical (some analytic) rj/s 'upper bounds' for K = 2, 3, 4 
(red, green and blue respectively). 



To conclude, we hope our work stimulates further research into holographic studies of 
the rj/s ratio in particular, and the dynamics of Lovelock theories of gravity and holography 
in general. It might be useful to explore further the effective potentials. These relatively 
simple functions seem to encode a great deal of information about the hypothetic dual 
CFT, information which we are only beginning to extract. 
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A. Three-point function parameters 

Conformal symmetry is powerful enough to determine the form of the 3-point function for 
the stress-energy tensor up to five constants and this are further constrained by conservation 
laws allowing us to reduce the number of independent parameters to three, A, B and C. 
Further one finds that Ward identities relate the two- and three-point functions and so Ct 
can be expressed in terms of these three constants. 

n n d _ 2 (d-2)(d+l)A-2B-4dC 

° T = ~2 (d-l)(d+l) ( } 

where is the area of the unit {d — 2)-sphere. Also as the energy flux one-point function 
is just a quocient of a three- and a two-point function t 2 and t± should also be writable in 
terms of the three parameters A, B and C. This was done in [16] yielding 

2d {d -3){d + l)dA + 3(d - 1) 2 B - 4(rf - l)(2d - 1)C 
t2 ~ d-l {d-2){d + l)A-2B-MC ' 

d (d+ l)(2{d - l) 2 - 3(d - 1) - 3)A + 2(d - l) 2 (d + l)B - 4{d - l)d{d + 1)C 
U ~ ~d-l~ (d-2)(d+l)A-2B-MC ' 

Then the scalar, vector and tensor constraints coming from from the positivity of energy 
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can also be translated into 

d (d-3){d+l)A + 2(d- l)B-4(d- 1)C 
~d- 1 (d-2)(d+l).4-2£-4de 



> 



( d _ 3 )(d + + (3d - 5)B - 8(d - 1)C 
d (d-2)(d+l).A-2B-4de " ' 1 > 

B — 2C 

- Ma ("- i! » ( (l -2)(d + iM-2B-4de so - 

constraints that can be identified with those coming from causality in the boundary theory 
in order to obtain expressions for A, B, C. However in order to unambiguously determine 
the normalization of these coefficients we can just obtain them in terms of Cr,i2>i4 and 
plug the expressions for these into them. We obtain then 

A (^~ 1 ) 3 T [d] 1 ( 2dAV'[A] 

(d - 2)3 tt^- 1 (-A)*/ 2 V (d - 4) 2 + 1 J 

- [ ^i^m + ( , _ ilf + M _ , r|A , 

(d-1) 2 ry 1 /(d 3 -3d 2 + 3d-4) AT"[A] 1 , NI w/rA A 
' ! + -(2(d- 3)d + 3)T'[A] . 



2(d - 2) 3 vr^ 1 (-A) d / 2 ^ (d-4) 2 2 

B. A curve through the parameter space of Lovelock gravities 

Exploring the full parameter space for a generic Lovelock theory is fairly untractable. As 
such, we choose to consider a particular curve through parameter space, simple enough 
that we can determine where along it a given theory preserves causality and stability. This 
curve corresponds to such a choice of parameters that the defining polynomial of the A-th 
order theory becomes (we set L = 1 here for simplicity) 

T[ 5 ] = l+ 5 + A 5 2 + ...= (l-|) (l-£)* 1 , (B.l) 

where there is just one free parameter as 

~ (K - 1)A ,„ . 

= - iTa~ ' (B - 2 > 

to ensure that the coefficient of order g on the polynomial (the Einstein-Hilbert term) is 
actually one. All the Lovelock coefficients can be written in terms of this parameter, 

(1 + A)((A-2)A-A) 
2(A-1)A 2 

(1 + A)2(A-2)((A-3)A-2A) 
~~ 2{K - 1)2 A 3 ' (R3) 
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We call this from now on the maximally degenerated trajectory (MDT). In figure 11 we 
show the projection of our curve on the X-fi-plane for the particular case K = 5. For the 
range A G [— K, 0], the cosmological constant, corresponding to the root with minimum 
absolute value of T = 0, is given by A, and beyond that by A. The curve parameterized by 
A connects Einstein- Hilbert gravity (A = —1) with the maximally degenerated Lovelock 
theory (A = A = —K), henceforth denoted by the maximally degenerate point (MDP). 
Given this curve, we want to see where instabilities or causality violation may occur. We 
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Figure 11: Projection of the points corresponding to (B.l) in the A-^-plane for K = 5. The dashed 
line correspond to the projection of the (K — l)th order trajectory, and in blue the corresponding 
interval of A values for each part of the curve is indicated. 

will restrict ourselves to the range A £ (— K, —1) section for any value K, i.e., to the part 
of the trajectory connecting the Einstein-Hilbert point with the MDP. We are looking for 
solutions of c?(r) equals zero or one for the stability or the causality analysis respectively. 
These constraints can be written as polynomial equations in g, and the appearance of new 
roots is signaled by a change of sign of the corresponding discriminant. An important point 
worth noticing here is that we are just interested in roots, g+, in the interval of allowed 
values for the function g 

<7*e[A,0], (B.4) 

where A is the closest to zero root of T[A] =0. The zeros of the discriminant separate the 
theory space spanned by (d, K) into distinct regions. By examining the properties of curve 
for a particular pair (d, K) we are guaranteed to have found the generic behaviour of the 
curve within the region containing that point. 

B.l Stability analysis 

All three potentials are bound to yield zero at the horizon and one at the boundary. Thus, 
any root of cf[g] = coming from outside the interval (B.4) has to enter it through the 
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horizon (g = 0) and in that case the horizon (which is always a zero for the equation) 
becomes a degenerate root. This is exactly what we studied when analyzing stability at 
the horizon - we look for the vanishing of the first coefficient in the horizon expansion. 
The same phenomenon occurs at the boundary for the causality analysis. Conversely, if 
any root leaves the interval between the boundary and the horizon one will also see a 
degenerate root appearing. The only other option for an instability to show up is as a new 
double root appearing in the bulk. In any case, it is clear that a study of the discriminants 
of the polynomial equations tells us exactly when these double roots occur. By carefully 
studying the regions in parameter space where the discriminant changes sign, one can learn 
where the stable and/or causal regions lie in parameter space. While this program is very 
involved in the generic case, the particular trajectory in parameter space we have chosen 
is simple enough so that this analysis of stability and causality can be done for all three 
potentials. 

B.l.l Stability at the horizon 

The horizon stability analysis can be performed by combining the expressions (B.3) for the 
Lovelock couplings along the MDT, with the stability constraints found in (4.17)-(4.18). 
For each constraint we get a quartic equation on A and we can study the number of 
real roots by examining the corresponding discriminants. For the lower boundary of the 
stability wedge, corresponding to the tensor channel, we get 

A$ = -16(d - lf{K - l) 4 K 2 (d - 3K - l)(d - 2K - 1) x 

[(d(162d - 1529) + 3768)K 4 - {d(3d(A5d - 391) + 1796) + 4044)K 3 (B.5) 
+ (d(d(3d(9d - 56) - 530) + 3576) - 504)K 2 

-4(d - l)(d(d(2d - 21) - 16) + 300)K + 8(d - l) 2 ((d - 2)d - 12)] , 

where s stands for stability and the number to the corresponding helicity of the channel 
considered. The h indicates that we are just considering stability at the horizon, for the 
moment. 

For general values of K this discriminant vanishes for d = 2K + 1 and d = 3K + 1. 
For d < 3K + 1 there is an interval of values for A around A = — K for which the tensor 
potential becomes unstable. This is simply due to the fact that the end point of the MDT 
is outside the stability wedge for d < 3K + 1 . Then, by considering high enough d for fixed 
K, the tensor channel becomes stable along our curve, close to the MDP. 

We can proceed in a similar manner with the upper boundary of the naive stability 
region, given by the sound channel. The discriminant in this case is 

= -432(d - 2){d - lf{d - 2K - l)(d - K - \){d - K)(K - l) A K 2 x (B.6) 
{{{K - 8)K + 8)d 2 - 2(K((K - 13)K + 6) + 8)d + K(20 - K(23K + 2)) + 8) 
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and it yields zero whenever d = 2K + 1 and 



d± 



K 3 - 13K 2 + 6K± V / (K-2) 2 K 2 (K 2 + K - 1) + 8 
K 2 — 8K + 8 



(B.7) 



The relevant root for the part of the curve we are interested in is actually d- . For d < d- 
there is a unstable interval of A values strictly inside A £ (—A, —1) as can be seen in 
figure 12. For d = 2K + 1 the MDP is at one of the edges of this interval. In the plot we 
can also check that for d > d_ this unstable interval disappears leaving a perfectly stable 
trajectory until it reaches the MDP. The situation changes dramatically as we increase the 
order of the Lovelock theory though. This is true for K < 7 - for K > 7 the denominator 
of d- changes sign and this critical value becomes negative and thus irrelevant. In other 
words, for K > 7 the unstable interval remains for every dimension. These statements are 
illustrated by figures 12 and 13. 



Figure 12: K = 3 trajectory (in black) and d = 7,8,9,10,11,12,15,20 helicity stability con- 
straint at the horizon (in blue). The unstable region is above the blue line. All the roots are in the 
figure in this case. 

B.1.2 Full stability 

For K > 7 it would appear that the sound channel is stable beyond the forbidden interval. 
This turns out not to be the case. As one goes beyond this interval the instability moves 
deeper in the bulk, so that there is always an instability in this channel. The full stability 
analysis is performed as follows. First notice that the two functions involved in the effective 
potentials, F'(x) and F"(x), have particularly simple expressions along the MDT: 




0.35 



F'{x{g)) 



(K -l)(g- A)(A + 1) (2A 2 + K (-A 2 



+ gk + k + g)) 



(B.8) 



{h 2 + K{-A 2 +gA + g)f 



F"(x(g)) 



2(K - l)(g - A)A 2 (A + 1)(K + A) 2 ((A- - 1)A + g(A + 1)) 
(K 2 + K(-K 2 + gK + g)) A 



(B.9) 
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Figure 13: K = 8 trajectory (in black) and d — 17, 20, 25, 30, 50, oo helicity stability constraint 
at the horizon (in blue and dashed black for d — > oo). The unstable region is above the blue line. 
All the roots are in the figure in this case. 



We have removed the explicit dependence on the derivatives of the function g(r) by using 
the polynomial equation (2.11). For instance one gets: 

In this way we can analyze the potentials just knowing the range of values g can take 
and avoiding solving the polynomial equation explicitly. In order to complete the stability 
analysis, now in the bulk of the geometry, we have to check if there is any value of g £ [A, 0] 
for which any of the potentials takes negative values. The appearance of any instability 
will be signaled with a change of sign of the discriminant of the corresponding polynomial 
equation c? = 19 . This is where a minimum of the potential touches the zero axis yielding 
a new double root of the equation. It is then enough to look for zeros of the discriminant 
and check if there is instabilities in the different disconnected regions. 

The full discriminant in the helicity zero case can be nicely written in terms of the 
corresponding discriminant of the horizon condition, 



(d- 1) 12 (A + 1) 12 ' 



A.o = im?7a i\t?A so , (B.ll) 



and then the appearance of instabilities in the bulk is controlled by the same discriminant 
analyzed in the precedent section. The A dependent factor vanishes just at the maximally 
symmetric point as A = is never possible in the present case. From the previous analysis 
we know that there are two naively stable regions for d < d- and just a full naively stable 
trajectory as we go to higher dimensions. We have to consider separately the allowed values 
for A on both sides of the unstable interval. 



19 This equation is not really polynomial but rational. However the denominator of the expression doesn't 
play any relevant role in this discussion and so we will restrict ourselves to the polynomial equation given 
by the numerator 
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For K < 7 and any value of A in between the forbidden interval and the MDP there is 
an interval of values for g where the potential goes below zero. This instability disappears 
for d > d-, at the same dimension as the naive unstable interval does (see left figure 14). 
Notice there's no instability either on the other allowed region as can be seen on the right 
hand side of figure 14. For higher values of K the instability in the allowed interval for A 




Figure 14: Helicity zero potential for K = 3 and d — 7,8, 9, 10, 11, 12 to the left and to the right 
of the forbidden interval on A, A = —2.5 and A = —1.15 respectively. 

close to the MDP never disappears. No matter how high the dimensionality some values of 
g lead to negative values of the sound potential (see figure 15). As one enters the stability 
wedge the instability simply moves away from the horizon and out towards the boundary 
of the geometry. 




Figure 15: Helicity zero potential for K = 8 and d 
right of the forbidden interval on A, A = —2.5 and A 



17, 20, 25, 30, 50, 100 to the left and to the 
— 1.15 respectively. 



The same analysis can be done in the helicity two channel. There the full discriminant 
turns out to be 

A 12 (K + A) 12 



(/<) 



(B.12) 



s2 ~ {d-l) 12 (A + l) 12 ^ s2 

In this case we just have one naively stable region for A G (— K, —1) and the tensor potential 
has been found to be stable everywhere there. Then the tensor channel is free of instabilities 
as long as it is so at the horizon (see figure 16). 

Even though we haven't made any reference to it yet, we didn't forget the shear 
potential. The reason for this disregard is that the corresponding stability constraint is 



- 36 - 




Figure 16: Helicity two potential for K = 8 and d = 17, 20, 25, 30, 50, 100 and cosmological 
constant A = —4. We find causality violation, but no instabilities. 



always less constraining than the helicity zero and two ones. In fact, it amounts just to 
A < A c , or equivalently rj/s > 0, when the other two constraints are respected (see eq. 4.9). 

All the features described in this section for the maximally degenerated trajectory can 
be easily observed in the cubic Lovelock case as we increase the dimensionality (see figures 
17). 

B.2 Causality analysis 

Let's now try doing the same kind of analysis with the causality constraints. These are of 
the form 

T'[A] + rAY"[A] = (B.13) 

For the MDT trajectory, solving this equation yields 

K + 2(K-l)r 
1 + 2{K-\)t 

Depending on the channel r takes different values namely (see (4.11)), 

2(d- 1) 



T2 

n - (B.15) 



TO 



(d-3)(d-4) 
d- 1 
d-3 
2(d- 1) 



(d-3) 



In order to respect A* < we must have r > — 2 (K-i) or T < — Wk-T) ■ Then, the crossing 
point of the MDT with the helicity two constraint is between the Einstein-Hilbert point 
and the MDP, and the crossing with the other two constraints is before the Einstein- 
Hilbert point. The helicity two crossing happens at a particular value of the Gauss-Bonnet 
coupling, namely 

K-1 K(K + 2(K-1)t) 
A " 2K (A' + 2(^-1)75)2 ■ ( ' 
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Figure 17: All the different regions in the third order Lovelock case for d = 9, 11, 20 (and zoom), 
from top to bottom. In black, the maximally degenerated line considered in the 'full' causality and 
stability analysis. For d < 11 the region of this curve closest to the maximally symmetric point 
(black dot) is unstable whereas for d > 11 it becomes stable. 

For the helicity two constraint the coefficient r goes to zero as d increases, that is, the 
value of A at the crossing point (the maximal value of A regarding causality) is the value at 
the maximally symmetric point with an increasing extra factor that goes to one as d — > oo. 
Then the crossing point approaches the MDP as d increases. This is explained by the fact 
that the helicity two causality reduces to Ct > for d — > oo, and hence our trajectory can 
only cross it at the MDP itself. 

The full causality analysis is quite more complicated that the stability one and not 
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very enlightening. Here we will simply cite some partial results in the most relevant case, 
the helicity two one. In that case the equation c|(r) = 1 can be reduced to a polynomial 
equation on g with discriminant 

A c2 = A 20 (l + A) 12 (K + A) 12 (A - A,) A c2 . (B.17) 

The factor A C 2 is a complicated expression which has zeroes at points irrelevant for our 
analysis. In particular, when restricted to the part of the curve connecting the Einstein- 
Hilbert point with the maximally symmetric point (A G (— K, — 1)) the discriminant has 
just one relevant root at A = A*, exactly the point determined by demanding causality at 
the boundary. The situation is analogous for the other two helicities. Then, causality in 
the MDT reduces to causality at the boundary. In particular, close to the MDP the only 
constraint comes from the helicity two channel, as has been discussed. 

In the same way as for the stability analysis, all the previously referred features can 
be seen in the third order case (see figures 17). 
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